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I . INTRODUCTION 

^ During  the  past  decade,  an  increased  demand  for  long-range  naviga- 
tion aides  has  been  exhibited  by  commercial  and  military  interests. 
Fishing  fleets,  oil  exploration  companies,  and  global  airlines  are  among 
the  host  of  potential  users  for  which  very  accurate,  remotely  available 
navigation  information  is  imperative.  OMEGA  is  a long  range  radio  navi- 
gation system  that  has  shown  considerable  promise  in  such  applications 
and  in  fact  has  demonstrated  a world-wide  capability  with  accuracies  on 
the  order  of  1 to  2 nautical  miles. 

Although  OMEGA  has  gained  modest  acceptance  within  the  navigation 
community  (particularly  in  marine  applications) , it  is  a propagation 
limited  system  and  several  problems  persist.  Certain  changes  in  the 
daytime  and  nighttime  characteristics  of  the  transmission  medium  cause 
day  to  night  changes  in  the  received  navigation  signals,  inducing  what 
is  known  as  diurnal  shift.  In  addition  to  this  repeatable  (and  somewhat 
predictable)  phenomenon,  unusual  solar  activity  can  cause  severe  and 
rapid  changes  in  the  medium  which  also  affects  the  received  signal. 

Large  volumes  of  predicted  diurnal  shift  values  have  been  published 
which  reduce  errors  caused  by  diurnal  shift,  but  their  useage  is  cumber- 
some at  best  and,  of  course,  they  cannot  account  for  errors  induced  by 
high  solar  activity  or  other  anomalous  propagation  conditions. 

This  dissertation  presents  two  new  techniques  which  may  be  used  to 
compensate  for  propagation-induced  errors  in  the  OMEGA  navigation  system. 
Chapter  II  includes  a description  of  OMEGA  and  appropriate  navigation 
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schemes  and  is  entirely  tutorial  in  nature.  Certain  kinds  of  unusual 


propagation  conditions  and  the  data  collecting  scheme  for  this  study  are 
also  described.  Chapter  III  introduces  a methodology  which  enables  one 
to  compensate  for  diurnal  shift  effects.  This  technique  resulted  from 
the  combined  efforts  of  Dr.  R.  G.  Brown,  Iowa  State  University,  and  the 
author  over  the  past  year.  Chapter  IV  presents  an  additional  compensa- 
tion method  developed  entirely  by  the  author,  which  enables  the  OMEGA 
user  to  detect  and  compensate  for  anomalous  propagation  conditions  en- 
tirely on  an  on-line  basis.  Both  the  methods  of  Chapters  III  and  IV  are 
deduced  from  physical  interpretations  of  VLF  (Very  Low  Frequency)  wave 
propagation,  and  offer  improved  accuracies  in  current  OMEGA  techniques 
and  increased  insight  into  the  VLF  navigation  problem. 
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II.  THE  OMEGA  NAVIGATION  SYSTEM 

A.  System  Description 

OMEGA  is  a very- long- range,  VLF  radio  navigation  system  which 
transmits  at  several  assigned  frequencies  in  the  10  to  14  kHz  frequency 
band.  Global  coverage  is  provided  by  eight  transmitting  stations  located 
at  various  points  around  the  world.  Due  to  its  very  low  transmitting 
frequency,  OMEGA  provides  adequate  navigation  signals  at  longer  ianges 
than  other  systems  and  is  generally  insensitive  to  prevailing  weather 
conditions.  The  system  is  normally  used  in  a single  frequency  hyperbolic 
mode  although  three  frequencies  are  available  and  it  can  be  used  in  a 
ranging  mode  (see  Section  II. B). 

The  primary  navigational  signal  (10.2  kHz)  and  the  remaining  two 
signals  (11  -j  kHz  and  13.6  kHz)  are  transmitted  sequentially  as  coherent 
time-synchronized  continuous  wave  bursts.  Note  that  the  transmitted 
frequencies  are  in  the  exact  ratio  9:10:12.  Each  burst  is  of  about  one- 
second  duration  and  is  repeated  at  ten-second  intervals.  The  trans- 
mitting format  is  unambiguous  and  allows  each  station  to  be  identified 
by  a decoding  of  the  received  sequence  of  bursts.  Transmissions  are 
derived  from  a bank  of  cesium  frequency  standards  and  each  station  is 
controlled  as  a source  of  standard  signals.  An  expanded  discussion  of 
OMEGA  system  operation  can  be  found  in  (1)  - (4),  and  the  interested 
reader  is  referred  to  these  documents  for  additional  details. 
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B.  Navigation  Using  OMEGA 

A convenient  model  for  describing  VLF  radio  propagation  is  that  of 
wave  propagation  in  a waveguide  (4) . In  this  case  the  boundaries  of  the 
waveguide  are  the  lower  D region  of  the  ionosphere  and  the  surface  of 
the  earth.  One  of  the  immediate  consequences  of  such  a model  is  that  if 
one  resonance  "mode"  is  dominant,  then  phase  and  amplitude  should  vary 
regularly  as  a function  of  distance  from  the  transmitter  without  fluctua- 
tions due  to  interference  between  various  modes.  In  our  case,  this  con- 
dition occurs  at  frequencies  near  10  kHz  and  at  distances  greater  than 
1000  km  and,  in  fact,  measurements  have  shown  that  outside  a 1000  km 
"near-field"  zone,  each  transmitter  radiates  a stable  signal  pattern 
which  is  repeated  in  a radial  direction  at  approximately  30  km  segments 
from  the  station  (Figure  1).  These  repeated  segments  are  proportional 


Figure  1.  Phase  - distance  relationship  at  10.2  kHz. 
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to  distance  and,  therefore,  also  to  the  cumulative  phase  (or  phase  delay 
as  it  is  more  commonly  known)  between  the  transmitting  station  and  a 
receiver.  The  above  phase  delay/distance  relationship  plays  a key  role 
in  using  OMEGA  as  a navigational  tool,  since  the  phases  of  the  OMEGA 
signals  can  be  measured  using  VLF  receivers.  Of  course  distance  is  pro- 
portional to  the  total  phase  delay  while  the  user  is  only  able  to  measure 
phase  difference  modulo  2ir,  i.e.,  if  <J>c  = 2ir(n  + a)  is  the  cumulative 
phase  difference  (n  - integer),  then  2ira  is  the  phase  difference  measured 
by  the  receiver.  We  therefore  tacitly  assume  that  the  user  has  resolved 
any  lane  ambiguities  via  some  external  position  fix.  Thus  distance  is 
obtained  from  the  relation 

distance  = time  delay  • velocity  of  phase  propagation 


where  Tj  is  the  time  required  for  the  phase  front  at  some  frequency  fj 
to  travel  from  a transmitting  station  to  the  receiver,  v^  is  the  presumed 
known  phase  velocity  at  fj,  and  <J> ^ is  the  measured  phase  delay  in  cycles. 
Note  that  the  accuracy  in  obtaining  a measure  of  distance  is  limited  not 
only  by  the  ability  to  measure  the  signal  phase  but  also  by  the  validity 
of  the  presumed  known  phase  velocity  v^. 

It  was  mentioned  previously  that  OMECA  can  be  employed  in  either  a 
hyperbolic  or  ranging  mode.  In  the  ranging  mode  the  user  measures  the 


phase  of  a transmitted  OMEGA  signal  with  a receiver  which  is  equipped 
with  a very  accurate  internal  oscillator.  This  oscillator  has  been 
previously  synchronized  with  the  master  oscillators  at  the  OMEGA  sites 
and  the  difference  in  phase  between  it  and  the  received  OMEGA  signal 
represents  a measure  of  range  from  the  transmitting  station.  Range 
determinations  from  two  or  three  transmitters  are  used  to  compute  posi- 
tion. 

OMEGA  is  also  used  in  a hyperbolic  mode.  In  this  mode,  the  phases 
of  signals  from  pairs  of  stations  are  differenced  and  used  to  establish 
position  along  a line-of-position  (LOP).  LOP's  (shown  in  Figure  2)  are 
lines  of  constant  phase  difference  between  a signal  received  from  station 
A and  a signal  received  from  station  B,  and  therefore  also  are  a measure 


Figure  2.  Hyperbolic  geometry. 


of  constant  range  difference  from  stations  A and  B.  When  referenced  to 
the  surface  of  the  earth  these  lines  are  described  by  sets  of  spherical 
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hyperbolas  unique  to  each  pair  of  stations.  By  selecting  different 
pairs  of  stations  the  navigator  generates  several  LOP's,  the  intersec- 
tion of  which  is  his  position  on  the  earth.  To  facilitate  the  construc- 
tion of  LOP  lattices  necessary  for  navigation,  the  U.  S.  Defense  Mapping 
Agency  has  published  standard  charts  and  tables  for  this  purpose.  The 
popularity  of  the  hyperbolic  configuration  lies  in  the  simplicity  and 
low-cost  of  the  receiver  hardware.  Since  any  absolute  radio-frequency 
phase  shifts  are  common  to  all  the  OMEGA  stations,  it  follows  that  these 
shifts  are  removed  in  the  differencing  operation.  This  eliminates  the 
need  for  the  user  to  have  an  extremely  precise  reference  oscillator  and 
allows  a considerable  cost-saving. 

C.  Propagation  Uncertainties 

In  the  previous  section  we  noted  that  we  could  estimate  the  distance 
from  an  OMEGA  transmitter  to  a receiver  by  multiplying  the  measured  time 
delay  at  an  appropriate  frequency  (10.2  kHz  in  this  case)  by  an  assumed 
nominal  phase  velocity.  Obviously  as  the  true  velocity  varies  from  the 
nominal,  errors  in  position  occur.  One  error  of  this  type  is  what  is 
known  as  diurnal  shift . Diurnal  shift  is  caused  by  a change  in  the 
effective  reflecting  height  of  the  ionosphere  between  night  and  day, 
thereby  changing  the  ohase  velocity.  If  not  corrected  for  diurnal  shift 
causes  errors  on  the  order  of  10  nautical  miles  over  typical  propagation 
paths.  Fortunately  the  velocity  variations  are  somewhat  predictable  and 
large  sets  of  phase  correction  tables  (or  PPC's  as  they  are  known)  are 
available  for  compensation  of  varying  propagation  parameters.  The  PPC's 
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list  corrections  as  a function  of  user  location,  hour,  day,  month,  and 
selected  transmitting  station,  and  allow  navigation  with  accuracies  of 
about  1 to  2 nautical  miles. 

Along  with  the  logistics  problem  associated  with  the  use  of  the 
tables,  the  PPC’s  are  "open-loop"  corrections.  Any  anomalous  atmospheric 
conditions  which  affect  the  signal  phase  velocity,  cause  navigational 
errors  which  can  neither  be  detected  nor  corrected  by  use  of  the  tables. 
The  most  important  of  these.  Sudden  Ionospheric  Disturbances  (SID's)  and 
Polar  Cap  Absorptions  (PCA's),  mainly  occur  during  periods  of  high  solar 
activity  and  are  characterized  by  a lowering  of  the  reflecting  height  of 
the  ionosphere  (5) . The  resultant  increase  in  phase  velocity  can  cause 
navigational  errors  on  the  order  of  5-10  n.m.  Since  recovery  from 
these  events  is  slow  (tens  of  minutes  to  several  hours  for  SID's,  several 
days  for  PCA's),  single-frequency  OMEGA  performance  suffers  in  comparison 
with  other  navigation  schemes  during  these  unusual  conditions  (6). 

D.  Phase  Monitoring  and  Recording 

The  U.  S.  Coast  Guard's  OMEGA  Navigation  System  Operations  Detail 
(ONSOD)  maintains  monitoring  stations  throughout  the  world  and  collects 
phase  measurement  data  on  all  the  OMEGA  frequencies.  This  is  done  pri- 
marily to  verify  the  accuracy  of  the  PPC's  but  also  serves  as  a check  on 
the  operation  of  the  transmitting  stations.  The  phase  measurements  are 
recorded  via  plot  printers  and  then  read  at  the  rate  of  one  measurement 
per  hour  by  site  personnel.  The  entire  receiver/recorder  system  has  an 
rms  accuracy  of  roughly  2 centicycles  (cec) . 


9 


Strip-charts  from  two  monitoring  sites  (North  Dakota  - station  D, 
and  Hawaii  - station  C)  were  obtained  from  the  Coast  Guard  for  the  month 
cf  March  1975  and  comprise  the  bulk  of  the  data  base  used  in  the  present 
study.  Phase  measurements  (at  all  three  OMEGA  frequencies)  were  read 
from  the  charts  at  the  rate  of  one  per  hour,  corrected  for  gross  lane 
count,  and  transferred  to  punched  computer  cards.  Propagation  paths  for 
which  the  data  was  read  include  B-D  (Trinidad  - North  Dakota) , C-D 
(Hawaii  - North  Dakota) , D-C  (North  Dakota  - Hawaii) , A-D  (Norway  - North 
Dakota) , and  H-C  (Japan  - Hawaii) . The  notation  B-D  indicates  that  the 
signal  was  transmitted  at  station  B (Trinidad)  and  recorded  at  station  D 
(North  Dakota);  thus  the  recorded  phase  measurement  represents  the  total 
phase  delay  over  the  path  B-D.  The  data  used  for  each  path  is  10-31  March 
for  B-D  and  C-D,  10-22  March  for  D-C,  and  1-31  March  for  A-D. 

In  addition  to  the  2 cec  measurement  accuracy  limitation  mentioned 
earlier,  other  errors  appear  to  be  present  in  the  data.  For  the  March 
1975  data  analyzed,  errors  apparently  occur  in  the  C-D  data  due  to  a dead 
zone  in  the  recording  equipment  at  the  monitoring  site.  The  strip-chart 
paper  used  at  the  site  has  a recording  width  of  100  cec  (1  cycle).  Thus 
when  the  phase  of  a signal  "crosses  over"  from  one  cycle  to  the  next,  the 
pen  plotter  resets  itself  from  100  cec  to  0 cec  and  continues  on.  When 
this  happens  on  the  C-D  recording,  the  plotter  continues  to  track  the  100 
cec  value  for  some  time  before  switching  over  to  zero.  The  error  induced 
by  this  equipment  delay  ranges  from  .5  to  3 cec  and  was  not  corrected  for 
in  this  study.  This  is  mentioned  merely  to  point  out  that  errors  of  a 
certain  magnitude  may  be  expected  in  the  data  base.  It  is  not  a reflec- 


10 


tion  on  the  quality  of  the  recording  scheme  which,  after  all,  is  only 
designed  to  meet  a 2 cec  rms  accuracy  specification. 

The  paths  selected  in  this  study  afford  a variety  of  propagation 
conditions.  Path  B-D  is  predominantly  South  - North  while  D-C  and  C-D 
are  East  - West  and  West  - East  respectively.  Path  A-D  was  included 
since  it  traverses  the  Greenland  ice  cap  and  causes  some  unusual  propa- 
gation problems.  H-C  (Japan  - Hawaii)  was  originally  chosen  as  an 
additional  West  - East  path  but,  because  of  frequency  synchronization 
problems  at  the  transmitting  station,  it  was  not  used  in  the  final  anal- 
ysis. 


The  information  recorded  on  data  cards  included: 

a.  hour  of  day  in  Greenwich  Mean  Time  (GMT) 

b.  day 

c . month 

d.  path 

e.  phase  measurements  at  10.2,  1 1'^'»  an^  13.6  kHz. 

f.  strip  chart  error  associated  with  each  phase  measurement 

g.  "unusual"  data  flag 

h.  day /night/transition  flag 

The  day/night/transition  flag  indicated  whether  the  total  path  was  in 
light  (day) , darkness  (night) , or  a combination  of  both  (transition) . 
These  definitions  of  day  and  night  are  standard  and  are  used  in  this 
context:  throughout  this  paper. 
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III.  COMPENSATING  FOR  DIURNAL  SHIFT 

A.  Overview 

Diurnal  shift  in  the  measured  phase  of  the  OMEGA  transmissions  is 
the  most  common  cause  of  error  in  the  OMEGA  system.  The  use  of  PPC 
tables  effectively  minimizes  this  problem  but  the  procedure  is  awkward 
and  slow.  In  addition,  with  the  recent  advances  in  microprocessor  tech- 
nology, there  has  been  increased  interest  In  developing  on-line,  real- 
time OMEGA  navigation  schemes  which  do  not  depend  on  PPC  tables.  One 
approach  is  to  implement,  on-line,  a simplified  version  of  the  prediction 
algorithms  used  to  generate  the  published  tables.  Beavers  and  others  (7) 
have  recently  suggested  that  such  a technique  is  now  feasible.  The  other 
approach  is  to  compute  a fictitious  propagation  time  from  on-line  phase 
measurements  made  at  more  than  one  frequency,  and  Papousek  and  Reder  (8) 
and  Pierce  (9)  have  proposed  navigation  schemes  based  on  this  concept. 

The  present  chapter  expands  on  this  latter  approach  and  offers  an  algor- 
ithm for  computing  a propagation  time  that  is  based  on  all  three  phase 
measurements.  The  concepts  of  group  velocity  and  associated  group  delay 
play  an  important  role  in  the  development  that  follows,  so  a few  remarks 
about  wave  propagation  in  a waveguide  are  appropriate  at  this  point. 

B.  Relation  of  Group  and  Phase  Velocity 

For  simple  parallel-plane  waveguide  modes,  the  phase  ind  group 
velocities  vary  with  frequency  as  shown  in  Figure  3.  These  velocities 
have  the  physical  interpretation  as  the  velocities  at  which  the  wave's 
fine  structure  (phase)  and  envelope  (group)  propagate  (10).  The  height 
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of  the  guide  determines  the  cut-off  frequency,  so  the  height  will  also 

affect  the  phase  and  group  velocities.  For  example,  increasing  the 

waveguide  height  in  a simple  parallel-plane  waveguide  lowers  the  cutoff 

frequency  at  the  dominant  mode.  For  a wave  at  a frequency  above  cutoff, 

the  group  velocity  is  then  increased  and  the  phase  velocity  decreased. 

Furthermore,  for  simple  modes  the  geometric  mean  of  the  phase  and  group 

velocities,  v and  v , is  invariant  and  equal  to  the  free  velocity,  i.e., 
P 8 


/v~  v = 
P g 


V 


0 


FREQUENCY 

Figure  3.  Phase  and  group  velocities  for  parallel-plane  waveguide. 


To  get  a better  feel  for  the  physical  meaning  of  group  velocity 
consider  the  example  of  transmitting  information  with  a carrier, 

Eg  coswt , amplitude  modulated  by  a signal  whose  frequency  is  Af=Aw/2ir. 


f j r 


Such  a signal  would  be  represented  by 


E = Eq(1  + m cosA(ot)coswt 


where  u Is  the  angular  frequency  and  m the  modulation  factor.  This  can 
be  expanded  to  show  the  presence  of  the  carrier  and  side-band  frequencies 

(ID 

mE0 

E = EgCoswt  + — [cos (no  + Aw)t  + cos(w  - Aw)t] 


When  such  a signal  is  propagated  in  the  z direction  under  conditions 
where  the  phase  velocity  varies  with  frequency,  i.e.,  in  a dispersive 
medium,  the  resultant  wave  is  written  as 

mE0 

E = E^cos^t  - gz)  + — — {cos[(w  + Aio)t  - (6  + Ag)z] 

+ [cos(u)  - Au)t  - (E  - Ag)z]} 

where  g is  phase  lag  per  unit  length  and  z is  the  distance  traveled. 
Using  well-known  trigonometric  identities  E can  also  be  written  as 


E = Eq[1  + m cos(Awt  - Agz)]cos(cot  - gz) 


The  bracketed  part  of  the  expression  represents  the  envelope  of  the  wave 
and  we  see  that  the  envelope  travels  in  the  z direction  with  a velocity 

1 

I 

I 
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If  the  frequency  spread  is  small,  this  may  be  written  as  the  limit 


dm 

Vg  dB 


(2) 


The  idea  of  group  and  phase  velocities  is  sometimes  visualized  in 

terms  of  Figure  4.  A person  "running  alongside"  the  wave  packet  as  it 

progresses  through  space  would  see  the  envelope  slipping  backward  with 

respect  to  the  component  waves.  This  is  because  it  is  moving  forward  with 

the  group  velocity  v , which  is  less  than  the  phase  velocity  of  the  com- 

g 

ponent  waves.  If  6 is  plotted  as  a function  of  co,  the  phase  and  group 

velocity  may  be  interpreted  geometrically.  Figure  5 shows  such  a plot 

for  wave-guide  propagation.  We  see  that  the  slope  of  B/w  = 1/v^  is 

always  less  than  that  of  4^-  = 1/v  so  that  v is  always  greater  than  v 

du  g p g 

but  both  approach  the  free  velocity  of  the  guide,  v^,  as  to  approaches 


I 
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'FINE  STRUCTURE"  VELOCITY  - PHASE  VELOCITY  - 
ENVELOPE  VELOCITY  > GROUP  VELOCITY  » -jf- 

Figure  4.  Modulated  traveling  wave. 
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C.  Composite  and  Difference  OMEGA 


The  theoretical  existence  of  an  invariant  free  velocity  throughout 
the  OMEGA  frequency  range  led  J.  A.  Pierce  in  1968  to  suggest  a naviga- 
tional scheme  called  composite  OMEGA  (9) . Consider  a flat-earth  propa- 
gation model  as  shown  in  Figure  6.  At  the  OMEGA  frequencies,  the  phase 
and  group  velocities  are  nearly  equal,  one  being  slightly  greater  than 
the  free  velocity  and  the  other  slightly  less  (see  for  example,  (1),  (5), 


or  (12)). 


Thus  the  geometric  mean  is  approximately  equal  to  the  arith- 


metic mean,  i.e. , 
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T = ^ LL 

g Aw  u2  _ “l 


where  <J>^  and  are  the  total  phase  shifts  in  radians  over  the  propa- 
gation path  at  angular  frequencies  to  and  to  respectively,  and  T is 


the  group  delay  referred  to  the  mid-frequency 


<0j  + 0)2 


For  composite  OMEGA,  Pierce  used  as  the  group  delay 


T _ _3 H 

8 “ - fi 


where  f^  = 10.2  kHz  and  f^  = 13.6  kHz,  and  and  ^ are  the  total  phase 
shifts  in  cycles.  The  group  delay  is  then  referred  to  the  mid-frequency 


f3  + fl 


= 11.9  kHz. 


The  phase  delay  can  also  be  referred  to  11.9  kHz  through  the  approxima- 


Vf3  + Vfl 


where  T and  T represent  the  phase  delays  at  f.  and  f,  respectively. 

P3  pl  31 

Substituting  Equations  4 and  5 into  Equation  3 then  yields  (after  some 
algebraic  manipulations),  Pierce's  composite  time  which  should  be  rela- 
tively insensitive  to  changes  in  the  ionospheric  height.  Note  that  T^  is 
a linear  combination  of  the  phase  delay  times  at  the  upper  and  lower 
OMEGA  frequencies. 
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T 

c 


*1 

« -jr-  + 2.25 

*-  1 


+ 2.25(T  - T ) 

P3  P1 


(6) 


The  way  in  which  the  composite  time  is  formed  can  also  be  viewed 
from  another  perspective.  In  the  discussion  above  we  developed  the  idea 
of  a composite  time  by  giving  equal  weight  to  the  phase  and  group  delays. 
If  instead  we  view  as  a weighted  sum  of  the  phase  and  group  delays 
referred  to  11.9  kHz,  i.e.  , 


T = kT  + (1  - k)T 
c g P 


it  works  out  that  T^  can  again  be  written  as  a linear  combination  of  the 
phase  delay  times  at  the  upper  and  lower  OMEGA  frequencies 


T = T + m(T  - T ) 


1 


P3  P1 


(7) 


irk(f3  + V l 

where  m=j|— ^ ^ ^ + 1 J . We  see  that  changing  the  "blending"  (k  factor) 
of  the  phase  and  group  delays  results  in  different  values  of  m and,  in 
fact,  both  extremes  in  blending  and  all  variations  in  between  can  be  ob- 
tained by  varying  m from  0.5  to  4.0  (Figure  7). 

The  concept  of  viewing  propagation  delay  as  a blending  of  group  and 
phase  delays  can  also  be  applied  to  other  navigation  schemes.  In  differ- 
ence OMEGA  the  propagation  delay  is  obtained  by  subtracting  the  phase 


delays  at  the  upper  and  lower  frequencies  and  dividing  by  the  difference 
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ALL  WEIGHT  ON 
PHASE  DELAY 
AT  11,9  KHz 


EQUAL  BLENDING 
OF  PHASE  AND 
GROUP  DELAYS 


ALL  WEIGHT  ON 
GROUP  DELAY 
AT  11.9  KHz 


2,25 


Figure  7.  Interpretation  of  Pierce's  m factor  in  terms  of 
blending  mid-frequency  phase  and  group  delays. 


frequency.  That  is, 


T _ 13 Zl 

d - f3  - ^ 


But  this  can  be  shown  to  be  precisely  the  composite  time  with  m=4,  that 
is  with  all  the  weight  on  the  group  delay.  So  the  difference  OMEGA 
delay  can  be  thought  of  as  a finite-difference  approximation  to  the  true 
group  delay  at  11.9  kHz.  In  the  OMEGA  case,  the  error  in  using  the 
approximation  is  less  than  .01  per  cent  (13).  Furthermore,  this  "error" 
is  not  one  of  substance,  it  simply  means  that  the  envelope  velocity  ob- 
tained from  the  finite-difference  form  corresponds  exactly  to  a group 
velocity  evaluated  at  a frequency  near  to,  but  not  exactly  at,  the  mid- 
frequency. This  is  assured  by  the  mean-value  theorem  of  calculus. 

In  summary,  we  see  that  difference  OMEGA  is  dependent  entirely  on 
group  delay  at  the  mid-frequency  of  11.9  kHz,  whereas  composite  OMEGA 


T 


employs  a blend  of  both  phase  and  group  delays  at  the  same  mid-frequency. 
Of  course,  the  mid-frequency  phase  and  group  delays  are  analytic  in  both 
cases,  because  they  are  formed  as  linear  combinations  of  measured  phase 
delays  at  the  end  frequencies.  Nevertheless,  this  mid-frequency  phase 
and  group  delay  viewpoint  is  useful  in  assessing  variations  in  delay 
times  due  to  changes  in  the  transmission  medium. 

The  result  of  implementing  composite  OMEGA  with  equal  weighting  of 
phase  and  group  delay  is  shown  in  Figure  8.  Propagation  time  is  plotted 
against  time  of  day  (GMT)  for  the  path  B-D  (Trinidad  to  North  Dakota) . 
Each  solid  line  represents  a time  span  of  one  day  so  the  entire  graph 
contains  22  days  of  superimposed  data. 

In  viewing  Figure  8,  we  see  that  a difference  in  the  day/night 
levels  of  about  25  gsec  is  quite  evident.  The  reason  for  this  diurnal 
shift  lies  in  the  fact  that  the  invariance  of  the  average  of  the  phase 
and  group  delays  is  only  valid  at  the  mean  height  of  energy  flow  in  the 
waveguide  (9).  In  a curved  waveguide  such  as  the  earth-ionosphere  guide, 
this  mean  height  is  located  somewhere  between  the  earth's  surface  and 
the  lower  ionosphere  (Figure  9).  Thus  from  purely  geometrical  consider- 
ations we  see  that,  for  some  fixed  distance  along  the  ground,  the  path 
at  a height  h'  increases  from  day  to  night.  Therefore  although  the  free 
velocity  is  invariant  at  h',  the  composite  time  as  seen  from  the  ground 
changes  from  day  to  night  because  the  mean  path  length  at  h'  changes. 
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D.  Three-Frequency  Difference  OMEGA 

Although  some  researchers  have  suggested  increasing  m to  compensate 
for  the  diurnal  shift  in  composite  OMEGA  (8),  ( 1 A ) , a reexamination  of 
Figure  9 suggests  another  approach.  Earlier  we  mentioned  that  group 
velocity  increases  as  we  go  from  the  day  to  the  night  portion  of  the 
waveguide.  One  might  speculate,  then,  whether  or  not  there  is  a fre- 
quency at  which  the  diurnal  increase  in  path  length  at  h'  is  exactly 
offset  by  a corresponding  increase  in  group  velocity,  thereby  resulting 
in  an  invariant  group  delay. 


Figure  9.  Curved-earth  model. 

Theory  indicates  that  such  a frequency  does  exist.  In  Figure  10, 

v /c  is  plotted  as  a function  of  frequency  for  several  values  of  iono- 
8 

spheric  height  (13).  We  see  that  the  curves  for  the  nominal  day  (h=70  km) 
and  nominal  night  (h=90  km)  conditions  cross  over  in  the  frequency  range 
12.0  to  12.5  kHz.  That  is,  the  group  velocity  during  the  day  and  the 
group  velocity  at  night  should  be  equal  at  some  particular  frequency. 

Not  only  is  the  existence  of  such  a frequency  of  interest  but  also 
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Figure  10.  Theoretical  group  velocities  for  different 
ionospheric  heights. 

whether  it  is  invariant  under  a variation  of  transmission  oarameters, 
e.g.,  direction  of  propagation,  ground  conductivity,  etc. 

To  compute  a group  delay  that  can  be  referred  to  any  frequency,  we 
utilize  all  three  OMEGA  frequencies.  Assume  we  have  three  time  delay 
measurements  T^ , T^,  T^  corresponding  to  the  three  OMEGA  frequencies  Wj 
^2 » ^2'  Each  of  these  time  delays  represents  a ratio  of  total  phase 
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shift  to  frequency,  i.e.. 


T 


i 


i 


2,  3 


It  is  tacitly  assumed  that  lane  ambiguities  for  all  three  frequencies 
have  been  resolved.  In  effect,  the  three  phase  measurements  give  us 
three  points  on  the  3d  versus  u>  curve  as  shown  in  Figure  11.  Now  assume 


Figure  11.  Phase  shift  versus  frequency. 


that  3d  can  be  approximated  as  a quadratic  function  of  w over  a reason- 
able range  of  u> ; i.e.,  let 


3d  = C2w2  + CjU  + CQ  (8) 

We  can  then  choose  the  coefficients  C^,  such  that  Bd  goes  through 
the  measured  4>^,  $3  points  as  shown  by  solving  the  equations 


I T- 


!l 
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C0  + C1W1  + Vl  = *1 


C0  + C1U2  + C2U2 


C0  + C1W3  + C2U3 


Returning  to  Equation  8 we  note  that 


T = d_  = d 
g v 


dto 


8 d6 


d(Bd) 

dm 


2C2w  + Cl 


(9) 


Substituting  the  solutions  for  and  C into  Equation  9 and  recognizing 


that  the  ratios  and  occur  in  the  exact  ratio  9:10:12,  the 


following  equation  results  (15), 


T = (60^—  - 66)T.  + (-100^—  + 105)To  + (40^  - 38)T- 
g 1 2 io„  3 


= KlTl  + K2T2  + K3T3 


(10) 


Equation  10  enables  one  to  compute  a group  time  delay  referred  to  any 
arbitrary  frequency  w,  in  terms  of  the  three  measured  phase  delays  T , 


T„,  and  T„.  Note  that  T is  a linear  function  of  T, , T„ , and  T„  and 
2 3 g 1 2 3 


that  the  sum  of  the  coefficients  (weight  factors)  is  unity. 

In  Equation  10  it  is  of  interest  to  look  at  variations  in  the  co- 
efficients of  Tj,  T2>  and  T3  with  frequency.  These  variations  are 
plotted  in  Figure  12  and  we  note  that  in  the  12.0-12.5  kHz  range  none  of 


the  coefficients  exceeds  6.  Purely  random  errors  in  the  phase  delay 
measurements  do,  of  course,  get  "amplified"  by  the  coefficients,  so 
large  values  are  undesirable.  How  undesirable,  though,  is  a matter  of 
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Figure  12.  Variation  of  K , K2 , coefficients  with  frequency. 

degree  and  depends  on  equipment  accuracy  and  the  application.  For 
example,  if  one  assumes  the  three  measurement  errors  associated  with  Tj , 
T2»  and  to  be  independent  and  each  having  unity  variance,  the  resul- 
tant rms  error  in  T would  be  as  shown  in  Figure  13.  It  should  be  ap- 
6 

parent  that  the  best  choice  of  reference  frequency  involves  a compromise 
between  the  noise  error  shown  in  Figure  13  and  the  diurnal  shift  error. 
Note  particularly  that  the  total  amplification  factor  ranges  from  5 to  8 
in  the  12.0  to  12.5  kHz  range. 
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Figure  13.  Normalized  rms  error  due  to  measurement  noise. 


An  example  of  applying  the  three-frequency  algorithm  to  one  path  is 
illustrated  in  Figure  14.  The  propagation  path  is  B-D;  thus  the  computed 
group  delays  represent  the  group  propagation  time  of  the  Trinidad  signal 
(station  B)  as  it  is  seen  at  North  Dakota  (station  D) . Each  of  the  four 
figures  in  Figure  14  contains  the  group  delays  for  22  days  computed  at  a 
particular  frequency.  Frequencies  of  11.7,  11.9,  12.2,  and  12.4  kHz 
were  selected  for  comparison.  In  the  graphs,  the  hours  0200-0900  corre- 
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spond  to  "night"  or  the  time  in  which  the  entire  propagation  path  is  in 
darkness;  and  1400-2100  corresponds  to  "day"  or  that  time  in  which  the 
entire  path  is  in  sunlight. 

In  viewing  Figure  14  several  items  are  of  interest. 

(1)  We  note  that  a diurnal  variation  is  quite  noticeable  at  11.7  kHz, 

decreases  somewhat  at  11.9  kHz,  and  is  essentially  zero  at  12.2  kHz. 
The  average  group  delays  during  the  day  and  night  periods  are 
plotted  as  a function  of  frequency  in  Figure  15,  and  we  see  that 
the  average  diurnal  shift  is  indeed  minimized  at  the  crossover  fre- 
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Figure  15.  Average  group  delay  versus  frequency.  Trinidad  - 
North  Dakota  (B-D) . 10-31  March  75. 


quency  of  12.2  kHz.  This  is  consistent  with  the  theoretical  curves 
in  Figure  10.  Assuming  a nighttime  height  of  90  km  and  a daytime 
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height  of  70  km,  the  group  velocity  at  frequencies  below  crossover 
is  higher  at  night  than  during  the  day;  therefore,  the  night  group- 
delay  times  are  lower  than  those  during  daytime. 

(2)  From  Figure  14  we  also  note  that  during  the  transition  periods 
between  night  and  day,  the  group  delay  decreases  sharply  and  then 
recovers.  Again  in  terms  of  Figure  10,  this  indicates  that  the 
ionospheric  height  is  at  some  level  between  70  and  90  km.  At  the 
crossover  frequency,  for  example,  the  group  velocity,  which  is  the 
same  for  night  and  day  moves  higher  ("towards"  80  km)  when  entering 
the  transition  period  and  then  returns  to  the  crossover  value. 

(3)  Finally,  we  see  that  the  noise-like  excursions  at  night  (random  or 
otherwise)  become  greater  as  the  frequency  increases  above  11.9  kHz. 
In  fact  an  examination  of  Figure  16  shows  that  the  standard  devia- 
tion of  the  group  delay  is  minimized  at  12.2  kHz  during  day  condi- 
tions and  at  11.8  kHz  at  night.  This  indicates  two  things: 

(i)  noise  amplification  which  is  of  little  consequence 
during  day  increases  in  importance  at  night,  and 
(ii)  phase  changes  due  to  small  fluctuations  of  ionosphere 

height  are  minimized  at  different  frequencies  during  day 
and  night.  This  is  to  be  expected  since  from  Figure  10 
we  note  that  the  90  (night)  and  80  km  curves  intersect 
at  a frequency  much  less  than  11.9  kHz  while  the  70  (day) 
and  80  km  curves  approach  each  other  as  frequency  in- 
creases. Thus  the  frequency  at  which  the  "noise"  in  the 
group  delay  is  minimized  is  a compromise  between  11.9  kHz, 
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where  the  instrument  error  amplification  is  least  (from 
Figure  13),  and  the  frequency  associated  with  minimizing 
random  fluctuations  in  the  transmission  medium. 


11.5  12.0  12.5 

frequency  (kHz) 

Figure  16.  Group  delay  standard  deviation  versus  frequency. 

Trinidad  - North  Dakota  (B-D) . 10-31  March  75. 

For  comparison  purposes,  propagation  times  computed  from  Pierce's 
composite  OMEGA  algorithm  are  shown  for  various  m factors  in  Figure  17. 
Again  propagation  time  is  plotted  against  time  of  day  (GMT)  and  it  is 
quite  apparent  why  increasing  the  m factor  decreases  the  diurnal  shift. 
For  Pierce's  original  value  of  m=2.25,  there  is  very  little  noise  asso- 
ciated with  the  composite  delay,  but  a diurnal  shift  of  about  30  usee  is 


Composite  time  versus  hour  of  day  (GMT)  for  various  m coefficients.  Trinidad 
North  Dakota  (B-D) . 10-31  March  1975.  Note:  The  figure  labeled  "Group 

Delay  - 11.9  kHz"  is  the  same  as  "Composite  Time,  m=4." 


quite  evident.  As  the  value  of  m is  increased  from  2.25  to  4.0  (that 
is,  as  more  weight  is  given  to  group  delay  and  less  to  phase  delay  when 
forming  the  composite  time),  the  average  diurnal  shift  reverses  itself. 
The  result  is,  for  this  path,  the  diurnal  shift  is  nullified  at  m=3.3. 

An  increase  in  noise  level  is  also  seen  as  m increases.  This  is  to  be 
expected  because  m is  a direct  noise  amplification  factor,  just  as  are 
the  Kj , K^,  coefficients  in  the  three-measurement  group  delay  formula. 

E.  Application  to  Hyperbolic  Navigation 

In  the  three-frequency  difference  OMEGA  development  we  tacitly 
assumed  that  we  were  using  OMEGA  in  a ranging  mode.  That  is,  the  group 
delays  were  computed  as  a linear  combination  of  phase  measurements  of 
signals  received  from  one  transmitting  site  (say,  for  example,  station 
A) . Then 


Vl, 


V2. 


K3T3 


where  are  constants  and  » T2  » T^  are  the  measured  phase 

AAA 

delays  from  station  A.  But  the  above  algorithm  can  also  be  applied  in 
a hyperbolic  situation  where  we  have  access  to  the  differences  in  phase 
between  signals  transmitted  by  a pair  of  stations.  Then,  applying  the 
three-frequency  scheme  to  the  phase  differences  between  stations  A and  B. 


T = K (T  - T,  ) + K.(T  - T„  ) + K,(T.  - T.  ) 

gA-B  1 LA  2 2 A 2B  3 3A  3B 


(K  T + K T + K T ) - (K  T + K T + K T ) 
A A J JA  B B B 


= T - T 
8A  8b 


So  the  group  delay  computed  with  phase  differences  from  stations 
A and  B is  equivalent  to  the  difference  between  group  delays  from  A and 
B and,  thus,  can  be  used  for  hyperbolic  navigation.  Note  that  the  propa- 
gation delay  being  a linear  combination  of  the  phase  delays  insures  that 
the  algorithm  can  be  applied  in  a hyperbolic  as  well  as  ranging  mode. 


F.  Summary  of  Results  for  the  Three-Frequency  Algorithm 


A complete  set  of  computer  generated  plots  which  reflect  the  re- 
sults of  applying  the  three-frequency  algorithm  to  the  test  data  is  given 
in  Appendix  A.  The  results  of  implementing  other  OMEGA  algorithms  are 
also  included  for  comparison.  Propagation  paths  included  are  Trinidad  - 
North  Dakota  (B-D) , North  Dakota  - Hawaii  (D-C) , and  Hawaii  - North 
Dakota  (C-D) . The  algorithms  presented  include  three-frequency  differ- 
ence OMEGA  referred  to  selected  frequencies,  composite  OMEGA  for  three 
values  of  m,  single-frequency  OMEGA  with  PPC  table  corrections,  and  a 
"modified"  composite  technique,  which  utilizes  nonlinear  combinations  of 
the  phase  measurements  (8).  Note  that  three-frequency  difference  OMEGA 
with  f = 11.9  kHz  is  equivalent  to  what  is  known  as  "difference  OMEGA"  and 
also  to  the  composite  OMEGA  algorithm  with  m-4.0.  Also  note  that  path 
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Chapter  II).  Non-reciprocal  path  effects  caused  by  the  earth's  magnetic 
field  are  also  in  evidence.  "Nominal"  propagation  times  are  plotted 
using  the  average  B-D  values  for  reference  and  we  see  that  the  D-C  group 
delays  are  a little  low  while  those  for  C-D  are  somewhat  high.  This  can 
be  predicted  from  theory  (12)  in  that  the  phase  velocities  for  East-West 
paths  (D-C)  are  slightly  higher,  and  for  West-East  paths  (C-D)  slightly 
lower  than  that  for  North-South  paths  (B-D) . 

In  terms  of  selecting  an  appropriate  reference  frequency  let  us 
examine  Figure  19.  For  each  of  the  three  paths  under  investigation  the 
rms  error  is  minimized  in  the  region  12.0-12.2  kHz.  Factors  which  in- 
fluence the  minimum-error  frequency  include  diurnal  shift,  instrument- 
error  amplification,  and  direction  of  propagation.  In  the  present  case 
the  amplification  factor  plays  the  dominant  role.  For  the  D-C  and  C-D 
paths,  it  interacts  with  the  directional  considerations  to  yield  a "best" 
frequency  of  12.1-12.2  kHz.  For  the  B-D  path  there  are  no  directional 
type  errors  and  the  trade-off  between  minimizing  noise  amplification  and 
diurnal  shift  results  in  a "best"  frequency  of  12.0  kHz. 

A summary  of  rms  errors  for  various  OMEGA  algorithms  is  given  in 
Table  I.  In  addition  to  the  difference  OMEGA  and  composite  algorithms 
previously  discussed,  a difference  OMEGA  scheme  compensated  with  PPC 
tables  (16)  is  included  for  completeness.  As  might  be  expected,  the 
algorithms  utilizing  PPC  corrections  were  the  most  accurate  followed  by 
composite  OMEGA  and  the  various  forms  of  three-frequency  difference 
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Figure  19.  RMS  error  versus  frequency  for  paths  B-D, 

D-C,  and  C-D.  10-31  March  75. 

OMEGA.  In  perspective.  Table  I illustrates  that  for  the  limited  data 
analyzed,  composite  and  difference  OMEGA  attain  similar  accuracies;  and 
while  neither  is  as  good  as  single-frequency  OMEGA  with  PPC  corrections 
the  errors  are  quite  acceptable  and,  hopefully,  not  as  subject  to  gross 
peak  errors  in  the  event  of  SID's. 


G.  Observations 


1.  It  has  been  shown  that  it  is  possible  to  compute  a group  delay 
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at  any  arbitrary  frequency.  Insensitivity  to  diurnal  phase  variations, 
random  noise  amplification,  and  other  factors  influence  the  selection  of 
a practical  operating  frequency. 

2.  For  each  propagation  path  examined,  a frequency  does  exist  for 
which  the  average  day  and  night  group  delays  are  equal.  Thus  it  appears 
that  diurnal  shift  can  be  minimized  by  a judicious  choice  of  reference 
frequency. 

3.  Navigational  accuracies  attained  using  three-frequency  differ- 
ence OMEGA  are  comparable  to  those  associated  with  composite  techniques. 

4.  The  phase/group  velocity  viewpoint  is  satisfying  from  a theoret- 
ical standpoint  and  offers  insight  in  developing  techniques  to  compensate 
for  other  propagation  phenomena,  e.g.,  SID's,  transition  effects,  etc. 


/ 
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IV.  DETECTING  AND  COMPENSATING  FOR  UNUSUAL  PROPAGATION  CONDITIONS 

A.  Overview 

One  of  the  major  limitations  of  OMEGA  is  that  it  is  susceptible  to 
anomalous  propagation  conditions  (17)  - (19).  SID/PCA  events  can  cause 
large  errors  in  single-frequency  OMEGA,  and  even  composite  techniques, 
which  have  reduced  sensitivity  to  SID's,  exhibit  degraded  accuracy  over 
unusual  paths,  e.g.,  the  Norway  - North  Dakota  (A-D)  path  which  traverses 
the  Greenland  ice  mass  (8) . In  this  chapter  a technique  is  presented  for 
compensating  computed  signal  propagation  times  by  sensing  changes  in  cer- 
tain transmission  parameters.  We  will  see  that  the  availability  of  phase 
measurements  at  three  frequencies  is  crucial  in  the  sense  of  providing 
information  about  propagation  conditions  along  the  path.  It  should  be 
noted  that  in  the  development  which  follows,  information  obtained  about 
propagation  conditions  can  be  utilized  in  any  navigation  scheme.  Thus, 
when  we  talk  about  propagation  delay,  it  is  with  respect  to  a delay  time 
referred  to  the  corresponding  velocity  (phase,  group,  or  otherwise)  with 
which  we  wish  to  navigate.  We  will  also  note  here  that  the  development 
assumes  that  the  user  is  operating  in  a direct  ranging  mode.  Extension 
to  hyperbolic  operation  depends  in  part  on  certain  geometric  constraints 
and  will  not  be  considered  below. 

B.  Using  Phase  Ratios  to  Sense  Propagation  Conditions 

To  begin  let  us  again  consider  the  navigation  scenario  where  we 
assume  that  we  are  in  a direct  ranging  mode.  Using  some  navigation 
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algorithm,  a user  computes  a propagation  delay  from  the  phase  measure- 
ments at  one  or  more  of  the  OMEGA  frequencies.  Then,  multiplying  the 
time  delay  by  some  assumed  nominal  velocity  of  propagation,  he  obtains 
a measure  of  distance  from  the  transmitting  station.  This  measure  of 
distance  deviates  from  the  actual  distance  depending  on  how  the  true 
velocity  differs  from  the  assumed  nominal  velocity.  That  is, 

d = T v 
comp  comp  nom 


and 


d = T v = T I — 
t comp  t comP\Vnom  ^ n0m 


so  the  propagation-induced  error  is 


error  = d^  - d 

t comp 


T V 

comp  nom  \ v / 
nom 


where  dt  = true  distance 


d = computed  distance 
comp 


v * true  velocity 


v * nominal  velocity  Cknown) 
nom  J 


T comp  = computed  time  delay 


(ID 


A 
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Note  that  the  error  is  a function  of  the  true  velocity  and  it  is  apparent 
that  any  method  which  allows  us  to  estimate  this  quantity  is  of  interest. 
For  example,  if  we  could  estimate  v during  a SID  period  then  presumably 
we  could  compensate  for  its  effects. 

To  see  that  the  ratio  of  phase  measurements  at  two  frequencies  can 
help  us  to  estimate  v , consider  the  case  of  propagation  in  a simple 
parallel-plane  waveguide  with  perfectly  conducting  boundaries.  The  ex- 
pression for  phase  velocity  in  this  case  is  (5) 


v = c{  1 
P 


(12) 


where  h is  the  ionospheric  height,  f the  frequency,  and  c the  velocity 
of  light.  In  Equation  12,  h is  the  only  unknown  parameter  but  we  cannot 
estimate  it  from  a phase  measurement  at  one  frequency  because  the  total 
phase  delay  also  depends  upon  distance,  an  additional  unknown;  i.e.. 


T,  = 


1 v x (h) 


(13) 


However,  if  a phase  measurement  at  another  frequency  were  also  available, 
then  presumably  we  could  "solve"  Equation  13  for  h and  then  obtain  an 
explicit  expression  for  v^.  Another  approach,  however,  is  to  look  at 
the  ratio  of  phase  measurements  at  two  frequencies,  i.e.. 


f 2 dVl  , . 

♦l  ' fl  T.  " £1  dv2  ' 


(U) 
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Distance  cancels  out  and  we  are  left  with  an  expression  in  which  h is 
the  only  variable.  Assuming  that  an  inverse  exists  we  can  write 


h = 


(15) 


or  equivalently,  since  v is  a function  of  h. 


(16) 


This  approach  is  useful  from  a practical  standpoint  in  that  the  precise 
functional  relationship  between  velocity  and  height  need  not  be  known. 

We  need  only  to  know  that  height  is  the  lone  variable  and  then  take 
sufficient  data  to  establish  an  empirical  model. 

Of  course,  in  the  case  of  OMEGA,  there  are  considerably  more  vari- 
ables at  play  than  just  ionospheric  height.  Height  is  the  primary  factor 
but  earth  conductivity,  direction  of  propagation,  ionospheric  conducti- 
vity gradient,  and  other  parameters  also  affect  the  propagation  velocity 
to  more  or  lesser  degrees.  Fortunately,  as  we  shall  see,  it  is  possible 
to  "lump"  many  of  these  parameters  into  generalized  variables  which  then 
can  be  viewed  as  functions  of  the  OMEGA  phase  ratios.  Let  us  first  con- 
sider the  approximate  model  for  phase  velocity  in  the  spherical  earth- 
ionosphere  waveguide  as  developed  by  Watt  (see  (5),  page  277).  Phase 
velocity  can  be  approximated  by 

Vp  = v^U  - ^'/a  + (1  - h/a)2[(2irn  + (J,^  + ^i)vQ/4iTfh] ^/2}  (17) 
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where  n = mode  number  (1  in  our  case) 

a = earth’s  radius 

h = reflecting  height  of  ionosphere 

h'  = mean  height  of  energy  flow  in  the  guide 

f = frequency 

Vq  = free  velocity  in  the  medium 
<J>  = ground  reflection  coefficient 


<J>^  = ionospheric  reflection  coefficient 
Thus  the  ratio  of  two  phase  measurements  can  be  written  as 


^l  fl  T1  fl  Vp2 


f.  {1  - h.'/a  + (1  - h/a)2[2ir  + (<b  + *.),v  /4Ttf  h]2/2} 

J, 1 S 2-i-Q * — (18) 

fl  - h^/a  + (1  - h/a)  [2it  + (<f>g  + 4>±) 2Vq/ 4-rrf 2h]  / 2} 


Over  the  OMEGA 
frequency  (5), 


frequency  range,  ^ and  are  essentially  independent  of 
so  Equation  18  can  be  written  as 


g(h,  hj , , (<{>g  + 4>1)) 


1 


(19) 


'77'"” 


i 
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It  also  turns  out  that  h'  can  be  viewed  as  a function  of  h and  (4>  + <f>  ) 

1 g 


so  that  Equation  19  can  be  rewritten  as 


IT  = g’(h,  K + 


(20) 


Equation  20  is  similar  to  Equation  14  except  that  we  now  have  an  addi- 
tional variable  cj>  +<)>..  With  the  OMEGA  transmissions,  however,  we  also 
g i 

have  at  our  disposal  another  independent  phase  ratio,  A similar 

kind  of  development  for  this  ratio  leads  finally  to  the  expressions 


if  4 R1  - 8"<h.  + *g> 


^ - R2  " g1 Chi  + *g> 


For  the  OMEGA  range  of  frequencies,  we  assume  g"  and  g'  are  invertible 
functions  and  we  thus  postulate  a model  for  estimating  velocity 
as 


v = f(RlS  R2) 


(21) 


where,  in  our  case,  R^  and  R2  are  measured  phase  ratios  and  the  model 
is  empirically  determined.  Note  that  in  the  development  above  (<f> . + <p  ) 
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C. 


Empirical  Determination  of  v as  a Function  of  R and  R 

8 i Z 


An  implementation  of  the  above  concepts  was  applied  using  three- 
frequency  difference  OMEGA  (reference  f requency=12 .0  kHz),  although  it 
should  be  noted  that  the  results  are  very  general  and  may  be  applied  to 
any  of  the  algorithms  discussed  in  Chapters  II  and  III. 

To  obtain  the  empirical  model  needed  to  estimate  the  true  group 
velocity  as  a function  of  R^  and  R^,  the  phase  measurements  from  paths 
B-D,  D-C,  and  C-D  were  used  as  raw  data.  The  phase  ratios  R^  and  R^ 
were  formed  from  each  measurement  and  since  the  distance  was  known  for 
each  path,  a measure  of  the  true  group  velocity  was  also  available. 

Thus  every  phase  measurement  was  represented  by  a point  in  three  dimen- 
sional space,  (v  , Rj , R^)  and  the  total  collection  of  these  three-tuples 
formed  the  data  base  for  the  model. 

Once  the  data  base  was  defined,  a fifth-order  polynomial  surface 
was  fit  to  the  points  using  standard  linear  regression  techniques  (20) . 
Group  velocity  was  designated  the  dependent  variable  and  an  expression 

was  generated  which  in  effect  estimated  v given  the  measured  values  of 

8 

Rj  and  R^,  i.e. , 


Vg  = a0  + “lRl  + “2R2  + a3RlR2  + *•••  + a20R25  (22) 

Values  for  coefficients  are  listed  with  their  respective  variables 
in  Table  II.  Those  variables  which  were  not  significant  at  level  .001 
were  not  included  in  the  basic  model  (Equation  22)  and  are  identified 
as  those  with  coefficients  equal  to  zero. 


Table  II.  Fifth-order  polynomial  fit  to  B-D,  C-D,  and  D-C  data 


Variable 

Coefficient 

Variable 

Coefficient 

R1 

-66.88 

R,V 

0.0 

R2 

-200.45 

Rl4 

-113.79 

R1R2 

0.0 

R24 

-3156.29 

Rl2 

0.0 

R1R24 

6365.45 

R22 

1360.59 

R14r2 

-90.61 

R!R22 

-2619.55 

„ 2„  3 

R1  R2 

0.0 

R12r2 

767.14 

riV 

0.0 

Rl3 

84.36 

Rl5 

0.0 

R23 

0.0 

R25 

0.0 

R1R23 

0.0 

Mean 

31.41 

r 3r 

0.0 

1 2 

The  surface  generated  by  the  regression  analysis  is  shown  in 
Figure  20.  Increasing  ionospheric  height  and  nominal  day  and  night  con- 
ditions are  indicated.  Interpreted  in  terms  of  the  group  velocity  ver- 
sus frequency  curves  in  Figure  10,  we  see  that  as  one  moves  along  the 
surface  from  nominal  day  (70  km)  through  transition  to  night  (90  km), 
the  group  velocity  starts  at  a fixed  level,  fluctuates  during  transition 
and  then  returns  to  the  fixed  level  at  night.  We  note  that  as  height 
decreases  below  70  km  (SID  disturbances)  the  surface  decreases  approxi- 
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This  is  because,  as  frequency  increases,  paths  "to  the  east"  (C-D)  have 

increasingly  lower  phase  velocities  than  paths  "to  the  west"  (D-C) ; 

thus  the  C-D  phase  ratios  are  higher  than  the  corresponding  D-C  ratios 

(1).  Note  that  the  estimated  v increases  as  we  move  from  D-C  to 

g 

C-D.  This  is  also  due  to  the  lower  phase  velocity  since  v^  is  inversely 

proportional  to  v^  and  therefore  increases  as  v^  decreases. 

The  surface  can  also  be  interpreted  in  terms  of  variations  in  the 

earth  reflection  coefficient  4>  . For  good  conducting  ground,  4>  is 

g g 

essentially  zero,  but  it  takes  on  negative  values  for  poorly  conducting 

ground  (5).  A decrease  in  <t>  results  in  higher  v (equation  17)  and  can 

g 8 

lead  to  group  velocities  which  occur  in  the  higher  region  "to  the  left" 
of  the  nominal  night  area.  Experimentally  this  is  confirmed  by  observing 
ratios  from  the  Norway  - North  Dakota  (A-D)  path  over  which  ground  con- 
ductivity is  extremely  low.  At  10  kHz,  in  fact,  ground  conductivity 
which  takes  on  values  from  3-5  mhos /meter  for  seawater  decreases  to 
something  on  the  order  of  2 x 10  ^ mho/meter  for  Greenland  ice  (1). 


D.  Experimental  Results 


In  this  section  the  effects  of  phase  ratio  compensation  are  dis- 
cussed as  they  concern  (i)  overall  system  accuracy  under  "nominal"  pro- 
pagation conditions,  (ii)  transition  periods,  (iii)  polar  conductivity 
paths,  and  (iv)  composite  applications.  With  the  exception  of  (iv) , 
the  three-frequency  difference  algorithm  was  selected  for  comparison  and 
the  following  procedure  was  used.  For  each  set  of  phase  measurements, 
the  phase  ratios  were  first  formed  and  an  estimate  of  group  velocity, 


v , was  obtained  from  the  fifth-order  model.  A propagation  delay, 
8est 

T , was  then  calculated  and  compensated  by  v /v  , i.e., 

8cal  8est  8nom 


8est  8cal\Vg 


where  v is  some  nominal  reference  velocity  and  T is  our  new 

8nom  8est 

estimate  of  propagation  delay  referred  to  v . For  this  study  v 

2 2 
nom  nom 

was  chosen  to  be  .9915  c (where  c is  the  velocity  of  light)  and  it 
appears  to  be  typical  of  day/night  velocities  for  the  B-D  path. 

(i)  Nominal  propagation  conditions.  The  results  of  applying  the 
compensation  to  the  three-frequency  difference  algorithm  (f = 12 . 0 kHz)  are 
shown  in  Figure  21  for  path  B-D.  Uncompensated  group  delays  are  plotted 
versus  time  of  day  (GMT)  in  the  first  figure  while  the  second  contains 
the  delays  compensated  as  discussed  above.  Random  variations  during  the 
day  and  night  periods  have  clearly  been  reduced.  The  severe  dips  charac- 
teristic of  the  transition  period  have  been  attenuated  but  not  eliminated. 
Such  "before"  and  "after"  pictures  for  all  the  paths  in\  atigated  (B-D, 
C-D,  D-C)  are  included  in  Appendix  B and  show  similar  results.  Table  III 
gives  the  overall  rms  accuracies  associated  with  both  uncompensated  and 
compensated  three-frequency  difference  OMEGA  and  we  see  that  there  is  an 
improvement  in  overall  accuracy  of  0.34-0.41  nautical  miles  depending 
upon  the  particular  path.  Comparing  Tables  I and  III,  we  note  that  with 
phase  ratio  compensation,  three-frequency  difference  OMEGA  achieves  a 


lower  rms  error  than  any  of  the  navigation  algorithms  which  do  not  use 
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PPC  tables.  Also  note  that  the  uncompensated  rms  errors  in  Table  III 
differ  slightly  from  those  in  Table  I because  of  different  reference 
velocities  used. 

(ii)  Transition  periods.  We  mentioned  earlier  that  the  transition 
times,  those  times  between  day  and  night,  appear  to  be  improved  very 
little  by  the  phase  ratio  correction  scheme.  In  fact,  by  deleting  a 
successively  larger  number  of  transition  hours  from  the  data  base,  we  can 
obtain  polynomial  fits  which  do  correspondingly  better  in  estimating  the 
group  velocity  from  and  R^.  This  can  be  seen  by  examining  the  square 
of  the  multiple  correlation  coefficient  associated  with  the  fit  analysis. 
This  number  essentially  indicates  the  amount  of  random  variation  in 

the  data  which  is  accounted  for  by  the  regression  model  (21) , and  we  note 
from  Table  IV  that  this  percentage  ranges  from  47%  in  the  case  where  all 
of  the  transition  times  are  considered  as  part  of  the  data  base,  to  72% 
when  the  transition  times  are  deleted.  This  indicates  that  transition 
effects  are  not  modeled  very  well  by  our  approach  and  suggests  that  other 
phenomena  (such  as  mode  energy  conversion)  are  at  work  during  this  period. 
We  should  note  in  passing  that  researchers  have  been  investigating  tran- 
sition phenomena  for  some  time  (22)  - (24) , and  transitional  effects  are 
still  not  well  understood  to  date.  Also  from  Table  IV  we  see  that  the 
day  and  night  delays  are  compensated  very  well  by  phase  ratio  compensation. 
In  fact,  with  the  compensation,  rms  error  for  path  B-D  is  reduced  to  0.6 
n.m.  during  day  and  night. 

(iii)  Polar  propagation  paths.  Another  interesting  possibility  for 
phase  ratio  compensation  is  the  application  to  propagation  paths  which 
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exhibit  very  low  surface  conductivity.  To  investigate  this,  a fifth- 
order  polynomial  was  generated  using  the  procedure  described  above  with 
path  A-D  (Norway  - North  Dakota)  included  in  the  data  base.  The  results 
for  the  low  conductivity  path  A-D  are  summarized  in  Table  V and  are  quite 
dramatic  in  that,  normally,  information  from  this  path  is  essentially 


Table  V.  Effect  of  phase  ratio  compensation  on  Norway  - North 
Dakota  (A-D)  rms  errors.  1-31  March  1975. 


Three-frequency  difference  OMEGA  (f-12.0  kHz) 
without  compensation 
with  phase  ratio  compensation 

Composite  OMEGA  (m=3.0) 

Modified  Composite  OMEGA 

Single  frequency  OMEGA 
(Using  PPC  tables) 

Difference  OMEGA 
(Using  PPC  tables) 


RMS  Error 


15.50  n.m. 
2.65  n.m. 

11.44  n.m. 

19.47  n.m. 


2.30  n.m. 


12.61  n.m. 


unusable  with  the  algorithms  which  do  not  use  PPC  tables.  Figure  22 
shows  the  results  of  using  this  particular  polynomial  to  compensate  for 
propagation  conditions  on  all  four  paths  (A-D,  B-D,  D-C,  and  C-D  respec- 
tively) . Although  performance  on  non-polar  paths  is  somewhat  degraded 
because  of  the  diurnal  shift  induced,  accuracies  on  all  four  paths  are 
within  2.7  nautical  miles  rms  which  is  quite  reasonable.  Note  that  the 
diurnal  error  occurs  totally  in  the  night  portion  of  the  guide  where  the 
phase  ratios  for  the  polar  and  the  "normal"  paths  tend  to  overlap.  This 
shift  in  the  night  values  occurs  because  the  estimated  Norway  - North 
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Group  delay  with  phase  ratio  compensation  versus  time  of  day  (GMT) 
for  various  propagation  paths.  10-31  March  1975. 
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Dakota  group  velocities  are  slightly  higher  than  for  the  other  paths, 
which  may  be  due  to  auroral  effects  (1),  or  mode  conversion  at  ice/ 
seawater  boundaries  (5) . 

(iv)  Composite  applications.  Phase  ratio  compensation  can  also  be 
applied  to  other  navigation  algorithms.  The  results  of  using  a fifth- 
order  polynomial  function  to  correct  composite  delays  (m=2.25)  are  shown 
in  Figure  23  (a)  and  (b).  Although  transition  effects  are  still  pro- 
nounced, overall  rms  error  has  been  improved  from  1.89  n.m.  for  (a)  to 
1 . 03  n.m.  in  (b) . 

E.  Application  to  SID  Events 

Sudden  Ionospheric  Disturbances  (SID's)  occur  in  conjunction  with 
high  solar  activity  and  are  characterized  by  sudden  phase  advances  in 
received  VLF  signals.  The  phase  advances  are  caused  by  an  increase  in 
signal  phase  velocity  due  to  a drop  in  the  effective  reflecting  height 
of  the  ionosphere  and,  since  they  are  not  predictable,  can  result  in 
severe  navigation  problems.  Since  SID's  can  be  thought  of  as  resulting 
from  ionospheric  height  changes,  the  phase  ratio  compensation  scheme 
presented  above  should  be  applicable  to  these  phenomena  also.  Unfor- 
tunately so  few  SID's  occurred  during  March  1975  that  utilization  of 
the  "usual"  data  base  was  not  possible. 

In  order  to  establish  some  sort  of  data  base  with  which  to  investi- 
gate this  application,  the  U.  S.  Coast  Guard  was  requested  to  furnish 
reduced  OMEGA  phase  measurements  from  the  Hawaii  and  North  Dakota  moni- 
toring sites.  This  data  (at  frequencies  of  10,2  and  13.6  kHz)  was  re- 
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ceived  and  appropriate  SID  events  identified  during  the  time  period 
1 January  1974  to  30  June  1975.  The  lack  of  11-^-  kHz  phase  data,  which 
was  not  furnished  because  of  certain  atmospheric  modeling  problems  which 
the  Coast  Guard  has  experienced,  required  a number  of  modifications  to 
be  instituted  for  the  analysis: 

(1)  It  was  decided  to  refer  the  group  delay  algorithm  to  11.9  kHz 
thereby  requiring  only  phase  measurements  at  two  frequencies. 

(2)  A linear  compensation  as  a function  of  one  ratio  was  assumed  in  the 
region  R^>1.3368. 

(3)  Each  of  the  three  paths  involved  (B-D,  C-D,  D-C)  were  examined  indi- 
vidually and  assigned  a different  value  for  the  slope  of  the  linear 
compensation. 

(A)  Only  SID's  which  occurred  during  "day"  periods  were  considered.  This 
restriction  essentially  forced  us  to  consider  only  those  values  of 
Rj  which  were  greater  than  about  1.3368,  i.e.,  in  the  region  where 
we  assumed  a linear  compensation. 

A complete  listing  of  the  days  and  paths  considered  for  the  analysis 
is  included  in  Appendix  C.  Appropriate  slopes  for  the  linear  compensa- 
tion functions  were  determined  from  an  examination  of  the  raw  phase 
ratios  and  are  also  tabulated. 

The  SID's  used  in  the  analysis  range  in  severity  from  very  minor  to 


potentially  catastrophic.  Results  are  summarized  in  Table  VI  and  we  see 
that  the  application  of  phase  ratio  compensation  does  improve  system 


60 


* 


61 


I 

I 

I 

I 


performance  even  with  the  rather  simple  algorithm  postulated  here. 


What's  more,  rms  figures  aside,  the  maximum  peak  errors  associated  with 


SID's  are  the  most  critical  insofar  as  navigation  is  concerned.  As  an 


example,  consider  the  effect  of  a particular  SID  as  shown  in  Figure  24. 
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Figure  24.  Propagation  error  versus  time  of  day  (GMT)  for  SID  event. 
Trinidad  - North  Dakota  (B-D) . 31  October  1974. 


The  SID  occurred  on  31  October  1974  along  the  propagation  path  B-D  and 


the  peak  error  incurred  by  the  single  frequency  OMEGA  algorithm  using 


PPC  tables  is  approximately  40  ysec.  or  about  7 nautical  miles.  Further- 


more, such  an  error  can  be  insiduous  in  the  sense  that  SID’s  cannot  be 
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predicted  and  the  user  may  not  be  aware  that  such  an  event  is  occurring. 
Also  from  Figure  24,  we  see  that  both  composite  OMEGA  and  the  compensated 
three-frequency  algorithm  do  creditable  jobs  in  mitigating  the  SID 
ef  fects . 

Another  possible  SID  application  of  the  phase  ratio  compensation 
technique  concerns  the  detection  of  the  SID's  as  they  occur.  In  Figure 
25,  smoothed  curves  of  propagation  errors  are  drawn  versus  the  ratio  , 
and  we  see  that  the  phase  ratio  can  be  viewed  as  an  indicator  of  whether 
or  not  a SID  is  in  progress.  That  is,  if  during  normal  operation  the 
ratio  crosses  a certain  threshold,  the  navigator  can  be  reasonably  sure 
that  a SID  is  occurring.  This  detection  capability  can  be  of  consid- 
erable value  in  many  applications  particularly  if  one  is  using  one  of 
the  PPC  table  algorithms  for  primary  navigation  (18). 

F.  Observations 

1.  Based  on  the  limited  amount  of  data  analyzed,  it  appears  that 
phase  measurement  ratios  can  be  used  to  sense  propagation  conditions  and 
correct  for  propagation-induced  errors.  System  performance  with  this 
form  of  compensation  is  comparable  to  that  of  single-frequency  OMEGA  com- 
pensated with  PPC  tables. 

2.  Results  appear  to  be  very  general  in  the  sense  that  the  compen- 
sation algorithm  can  be  made  independent  of  the  transmission  path. 

3.  Inherent  to  the  phase  ratio  technique  is  a capability  to  detect 
as  well  as  compensate  for  SID  events,  the  detection  aspect  of  which  is 
not  available  in  other  OMEGA  algorithms. 


V.  SUMMARY 


The  two  analytical  techniques  presented  in  this  dissertation  offer 
new  methods  for  compensating  for  propagation- induced  errors  in  the  OMEGA 
navigation  system.  Both  techniques  were  derived  primarily  from  consid- 
eration of  the  physical  phenomena  involved  with  VLF  propagation  and  are 
theoretically  satisfying  in  this  respect. 

The  three- frequency  difference  OMEGA  approach  illustrated  that 
errors  due  to  the  diurnal  phase  shift  of  the  OMEGA  signals  could  be  com- 
pensated for  by  a judicious  selection  of  the  frequency  to  which  the  group 
time  delay  was  referred.  Overall  navigation  accuracies  using  this  scheme 
were  comparable  to  those  of  composite  OMEGA. 

Another  technique,  phase  ratio  compensation,  was  used  to  sense  pro- 
pagation conditions  along  the  transmission  path  and  to  correct  for  pro- 
pagation-induced errors.  It  was  demonstrated  that  this  form  of 
compensation  significantly  improved  the  three-frequency  difference  OMEGA 
performance  and  in  the  case  of  anomalous  propagation  conditions  (i.e., 
SID's,  polar  transmission  paths)  allowed  navigation  when  one  or  more  of 
the  other  algorithms  was,  for  all  practical  purposes,  inoperative. 
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VIII.  APPENDIX  A:  PROPAGATION  TIMES  FOR  VARIOUS 

OMEGA  CONFIGURATIONS 
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Figure  60.  Uncompensated  and  compensated  group  delay  (f=12.0  kHz)  versus 
time  of  day  (GMT) . North  Dakota  - Hawaii  (D-C) . 10-22  March 
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Figure  61.  Uncompensated  and  compensated  group  delay  (£<*12.0  kHz)  versus 
time  of  day  (GMT).  Hawaii  - North  Dakota  (C-D) . 10-31  March 
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X.  APPENDIX  C:  DATA  BASE  FOR  SID  ANALYSIS 

Path  Day  Hour  (GMT) 

B-D  1 Feb  74  1900-2200 

B-D  3 Feb  74  1900 

B-D  4 Feb  74  2200 

B-D  5 Feb  74  2100 

B-D  15  Feb  74  1600,  1900 

B-D  16  Feb  74  1900-2200 

B-D  21  Feb  74  1600,  1700,  2200 

B-D  22  Feb  74  2100-2200 

B-D  13  Apr  74  1600-1900 

B-D  1 Jul  74  1900 

B-D  4 Jul  74  1500-2000 

B-D  5 Jul  74  1700,  2300 

B-D  6 Jul  74  1300-1500 

B-D  6 Jul  74  2000-2200 

B-D  8 Jul  74  1800 

B-D  2 Aug  74  1600-1900 

B-D  6 Aug  74  1800 

B-D  28  Aug  74  1200 

B-D  31  Aug  74  2000 

B-D  8 Oct  74  1500-2000 

B-D  10  Oct  74  1600-1900 

B-D  11  Oct  74  1800-1900 


L 


142 


i: 


i 

Path 

Day 

Hour  (GMT) 

• • 

f- 

B-D 

12 

Oct 

74 

1900 

I 

B-D 

18 

Oct 

74 

1300-1500 

1 

B-D 

30 

Oct 

74 

1700 

B-D 

31 

Oct 

74 

1600-1900 

u 

B-D 

3 

Nov 

74 

1800,  2200 

[ 

1. 

B-D 

5 

Nov 

74 

1600-1700 

B-D 

8 

Nov 

74 

1500-2100 

B-D 

11 

Feb 

75 

2000-2100 

B-D 

24 

Feb 

75 

1800 

LI 

D-C 

25 

Dec 

74 

2100-2200 

D-C 

5 

May 

75 

2200 

D-C 

19 

May 

75 

2100-2200 

D-C 

26 

May 

75 

2100-2400 

D-C 

27 

May 

75 

0100-0300 

r 

C-D 

10 

Mar 

75 

2200-2400 

C-D 

13 

Mar 

75 

2000-2400 

C-D 

14 

Mar 

75 

2000-2400 

p 

C-D 

15 

Mar 

75 

1900-2400 

C-D 

27 

Mar 

75 

2000-2100 

C-D 

9 

Apr 

75 

2000-2400 

I 

C-D 

21 

Apr 

75 

1900-2400 

1 

C-D 

22 

Apr 

75 

1900-2400 

